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Perturbations 



O 



o 

u 

o 



> 

m 

rn 
O 
m 



Taotao Qiu,i'2'El Xian Gao,3.4,5,[t| and Emmanuel N. Saridakisf^'^'l 

Leung Center for Cosmology and Particle Astrophysics National Taiwan University, Taipei 106, Taiwan 
''^ Department of Physics, National Taiwan University, Taipei 10617, Taiwan 
Astroparticule et Cosmologie (APC), UMR 7164-CNRS, Universite Denis Diderot-Paris 7, 
10 rue Alice Domon et Leonie Duquet, 75205 Pans, France 
Lahoratoire de Physique Theorique, Ecole Normale Superieure, 24 rue Lhomond, 75231 Pans, France 
^Institut d'Astrophysique de Pans (lAP), UMR 7095-CNRS, 
Universite Pierre et Mane Cune-Paris 6, 98bis Boulevard Arago, 75014 Pans, France 
^Physics Division, National Technical University of Athens, 15780 Zografou Campus, Athens, Greece 
^Instituto de Fisica, Pontificia Universidad de Catolica de Valparaiso, Casilla 4950, Valparaiso, Chile 

We investigate non-singular bounce realizations in the framework of ghost-free generalized Galileon 
cosmology, which furthermore can be free of the anisotropy problem. Considering an Ekpyrotic-like 
potential we can obtain a total Equation-of-State (EoS) larger than one in the contracting phase, 
which is necessary for the evolution to be stable against small anisotropic fluctuations. Since such 
a large EoS forbids the Galileon field to generate the desired form of perturbations, we addition- 
ally introduce the curvaton field which can in general produce the observed nearly scale-invariant 
spectrum. In particular, we provide approximate analytical and exact semi-analytical expressions 
under which the bouncing scenario is consistent with observations. Finally, the combined Galileon- 
curvaton system is free of the Big-Rip after the bounce. 

PACS numbers: 98.80.-k, 04.50.Kd, 98.80.Cq 



I. INTRODUCTION 

Non-singular bouncing cosmology [l[ has gained sig- 
nificant interest in recent studies of the early universe. 
The main reason for such a research direction is that the 
most popular paradigm of the early universe, namely in- 
flation, still suffers from the "Big-Bang singularity" prob- 
lem, which however can be naturally avoided in non- 
singular bouncing or cyclic cosmologies. Additionally, 
these paradigms can also solve the horizon, flatness and 
monopole problems, and make compatible observational 
predictions such as nearly scale- invariant power spectrum 
and moderate non-Gaussianities [l-ll]. Therefore, they 
are recently considered as good alternatives to inflation. 

In order to realize a successful bounce several require- 
ments must be fulfilled. First of all, the basic condition is 
to have the Hubble parameter change its sign from nega- 
tive to positive at the bounce, which implies that during 
the bouncing phase the Null Energy Condition (NEC) 
must be violated, with the total EoS of the universe go- 
ing below —1 0,0. The NEC violation in the context 
of General Relativity is nontrivial 0, usually leading to 
ghost degree(s) of freedom , which would demand ei- 
ther ghost-elimination mechanisms or an extended anal- 
ysis to a modified gravity context [lOt . 

Apart from the above basic condition, in order for a 
bounce to be a successful alternative to inflation it should 
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also solve the other Big-Bang problems, and moreover it 
should produce a nearly scale-invariant power spectrum 
as required by observations fllj . These impose more 
stringent constraints on the bounce evolution, especially 
in the contracting phase. For instance, the horizon prob- 
lem can be solved if the quantum fluctuations in the far 
past lie deep inside the horizon, while they should exit the 
horizon in the contracting phase in order to generate per- 
turbations compatible with observations, provided that 
inflation is absent in bouncing scenario. This requires a 
total EoS satisfying w > —1/3 in the contracting phase 
[isj . However, the scale-invariance of the perturba- 
tions is even harder to be achieved. In particular, as it 
was initially shown in Q , if the perturbations generated 
in the contracting phase are purely adiabatic, the EoS of 
the contracting universe should satisfy w « in order to 
produce the desired spectrum. 

However, although bounce models with total EoS w w 
before the bounce, namely the "matter bounce" , could 
lead to nearly scale-invariant power spectrum, they gen- 
erally suffer from the "anisotropy problem" [15j in the 
contracting phase. In particular, in 4D General Relativ- 
ity a tiny amount of anisotropic fluctuation from the sim- 
ple isotropic Friedmann- Robertson- Walker (FRW) geom- 
etry in the contracting phase, would increase as a~^, 
where a is the scale factor. Thus it would finally dom- 
inate over matter-like background, leading to a Big- 
Crunch singularity with complete anisotropy instead of 
a bounce, unless one impose a strong fine-tuning of the 
model parameters and the initial amount of anisotropy 
in order to obtain a bounce before the domination of the 
anisotropic term. In that sense the "matter bounce" sce- 
nario is not stable against cosmological anisotropy (for 
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its similar problem in the presence of radiation see jl6|). 
For this reason, we must construct scenarios with total 
EoS larger than 1 in order to prevent the dominance of 
anisotropy. However, as we mentioned above, a differ- 
ent EoS may not be able to provide the scale-invariant 
power spectrum, if we insist on applying the simple 
adiabatic mechanism of generating primordial perturba- 
tions^. Therefore, we should resort to alternative mech- 
anisms, such as adiabatic, entropy and conformal ones 

In the present work we investigate the bounce real- 
ization in the framework of recently proposed general- 
ized Galileon cosmology [12, (see also [13] for var- 
ious developments). Due to the delicate design of the 
Lagrangian form such a theory, which contains higher- 
order derivatives, can keep its equation of motion up to 
second-order and thus is free of ghosts (this was pioneered 
by the work by Horndeski [25| ) , but it can indeed provide 
extra degree(s) of freedom in order to violate NEC. Re- 
cently, in [2^ , the first ghost-free bounce model based on 
Galileon cosmology was constructed by one of the present 
authors and collaborators (see also [23), and hence in 
this article we will consummate this class of models by 
addressing the problems mentioned above. Note that al- 
ternative scenarios addressing the anisotropy problem in 
Galileon bouncing cosmologies have been presented in 
psj . of which before contracting with w > I, the uni- 
verse can be dominated by cold matter [2^ , where scale- 
invariant perturbations could be generated. 

First of all, by introducing an Ekpyrotic-like negative 
potential we can easily obtain a very large EoS in the 
contracting phase, thus the anisotropy problem will be 
eliminated. However, as mentioned above, a large EoS 
forbids the Galileon field to generate the desired form 
of perturbations, thus as a next step we additionally in- 
troduce the curvaton field which is suitably coupled to 
the Galileon field, such that the nearly scale-invariant 
spectrum can be produced. Finally, we perform a com- 
plete analysis of the behavior around the bounce point of 
the full Galileon-curvaton system, making use of the "in- 
verse" reconstruction procedure [30[ , showing that with a 
proper choice of the Lagrangian functional forms a non- 
singular bounce can be reconstructed, which can connect 
smoothly to the matter-domination era and moreover al- 
leviate the Big- Rip singularity which appears in j26j . 

The plan of the work is the following: In section [H] 
we briefly review the anisotropy problem. In section IIIII 
we present the bouncing background evolution before, 
during and after the bounce, and we show that the per- 
turbations are stable and free of ghosts. In section HVl we 



^ We would like to mention that here by "simple adiabatic mecha- 
nism" we mean that the perturbations generated are purely adi- 
abatic, and the EoS remains constant. However, the term "adia- 
batic mechanism" which was first proposed in [TtII in "Ekpyrotic" 
scenarios [l8l |. refers to the mechanism that generates adiabatic 
perturbations via varying EoS. 



analyze the curvaton mechanism that produces nearly 
scale-invariant perturbations. In section |V] we perform 
a semi-analytical procedure in order to reconstruct an 
exact bouncing solution that is not followed by a Big- 
Rip. Finally, in section I VII we summarize and we dis- 
cuss the obtained results. Throughout the manuscript 
we use the (— , metric signature, and units in 
which Mpi = I/VSttG = 1. 



II. THE ANISOTROPY PROBLEM 

The anisotropy problem is a notorious problem that 
generally exists in bouncing models with w < 1 in con- 
tracting phase [l^ . In General Relativity, if we allow the 
existence of a non-zero anisotropy at the beginning of the 
contraction it will evolve scaling as a~^{t). In order to 
demonstrate this more transparently, without loss of gen- 
erality we consider as an example the simple anisotropic 
Bianchi-IX metric [3]| : 

ds^ = -dt^ + a\t) J2 e^^'^*^dx'^ , (1) 

with /3i(t) 4- /32(t) + l33{t) = 0. The Friedmann Equation 
writes as: 

SH' = Pu + l(j2$f^ , (2) 

where H = a/a the Hubble parameter and with p„ in- 
corporating all the fluids in the universe. The /3i's satisfy 
the equations 

/3, + 3ff A = , (3) 

which provide the solutions pi oc a^^{t). Since the second 
part in the right hand side of equation ([2]) can be con- 
sidered as an effective anisotropy term , we conclude 
that 

x'^\(jZpi^^a-%t), (4) 

and thus the anisotropy term corresponds to an effective 
energy density with EoS w ~ 1. Although in an expand- 
ing universe this term is always sub-dominant and thus 
isotropization can be achieved, in a contracting case, as 
long as it is initially non-zero (even arbitrarily small), 
the anisotropy will grow fast and become dominant over 
all species with EoS less than 1, leading finally to a col- 
lapsing anisotropic universe. For this reason, in order 
to avoid the domination of a possible anisotropic fluctu- 
ation, one has to realize a contracting background that 
evolves even faster, which requires an EoS larger than 
unity in the contracting phase. 
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III. THE GALILEON BOUNCE 

In the previous section we briefly showed that in order 
to reahze a bounce we need an effective EoS w > 1 in the 
contracting phase, in order to avoid the domination of an 
anisotropic fluctuation. In this section we formulate the 
bounce realization in generalized Galileon cosmology. 

In the generalized Galileon scenario, where the coeffi- 
cients of the various action-terms are considered as func- 
tions of the scalar field, the corresponding action can be 
written as [2311: 



(5) 



i=2 



where 



£3 = -G3(<^,x)n0 

U^Gi{(^,X)R + GA,x[{U<j>f 
1 



- G^.x [{U4>f - m4>) (V^V,0) (V^V^c/)) 

+2(V^V„(^) (V"V^0) (V^V^0)] . (6) 

In this action the functions K and Gi [i = 3, 4, 5) de- 
pend on the scalar field 4) and its kinetic energy X = 
— ^V^^V^, while R is the Ricci scalar and G^v is the 
Einstein tensor. Moreover, Gi^x and Gi^^ {i = 3, 4, 5) 
denote the partial derivatives of Gi with respect to X 
and 0, (Gi_x = dGi/dX and G^.^ = dGi/dc/)), and the 
box operator is constructed from covariant derivatives: 
D(f> = g^'^V y4). In the following wc focus on the case 

if (0, X) = X- G3(0, X) = .g(0)X, 

G4(0,X) = i G5(0,X)=O. (7) 

Therefore, the action that we are going to use reads: 



S 



1 



-R 



1 
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(8) 



We now proceed to a detailed investigation of the above 
scenario. Firstly, in the following subsection wc provide 
approximate analytical solutions at the far past before 
the bounce, around the bouncing regime, and after the 
bounce, at the background level. Then in the next sub- 
section, we analyze the perturbation behavior. 



A. Background evolution: analytical results 

In the following we impose a flat Friedmann- 
Robertson- Walker (FRW) background metric of the form 



ds^ = —N'^{t)dt^ + a?{t)d'x^ , where t is the cosmic time, 
are the comoving spatial coordinates, N(t) is the lapse 

function, and a{t) is the scale factor. Varying the action 
with respect to N{t) and a{i) respectively, and setting 

A'^ = 1, we obtain the Friedmann equations 



H^ = \P^ H = -\{p + P). 



(9) 



Additionally, we have defined the effective energy density 
and pressure as: 



P =\^^- V{cp) - g4>^ , 
and thus the total EoS of the universe is just 

_P _ 4? - 2V{4>) - 25020 
^ ~ P ~ 02 4. 2^(0) + 6gi703 ' 



(10) 
(11) 

(12) 



Finally, variation of (|8]) with respect to the Galileon field 
provides its evolution equation: 



p0 + r0 + y^, = 



where 



V 

r 



1 



1'^ 



6gH<p 

igH4) ( 
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2H 



(13) 

(14) 
(15) 



In the following we will suitably choose the potential 
V{(j)) in order to obtain a very large positive equation of 
state in relation (|12p . namely w > 1, so that our model 
will not suffer from the anisotropy problem in the con- 
tracting phase. Meanwhile, when the Galileon term Gn0 
becomes more and more important, the EoS becomes 
negative and eventually triggers the bounce. 

As a specific example, we choose the potential to be 



(16) 



with Vb,c > 0; namely a negative exponential poten- 
tial, which is the usual one in Ekpyrotic scenarios (isj . 
Within this choice, when the nonlinear kinetic term is 
relatively small we obtain p < P and thus w > I. One 
could ask whether the negative potential would lead to a 
negative energy density, however the more negative the 
potential is, the steeper it is, and thus it gives rise to 
larger kinetic term, which can compensate the potential 
negativity (this will be verified later on). 

Let us proceed to a qualitative investigation of the dy- 
namics of such scenario. For simplicity we consider that 
increases from negative to positive during the evolu- 
tion, therefore at the beginning when starts at a large 
negative value the potential (|16|) is in its "slow-varying" 
region, and the field moves slowly. In this region the non- 
linear term Gn0 will have a small contribution to the 
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action. Similarly to the Ekpyrotic models, the universe 
will contract with a very large positive EoS, but as time 
passes (j) moves towards positive values and its velocity 
increases, the effect of the GDc^ term is enhanced, and it 
can trigger the bounce. However, after the bounce and 
due to the large slope of the potential, the field kinetic 
energy could increase to unacceptably large values which 
could spoil the validity of the effective theory. Therefore, 
we should need some mechanisms to slow the field mo- 
tion down, and this will be demonstrated in detail in the 
next sections. 

In the following we proceed to the quantitative investi- 
gation of the scenario, extracting approximate analytical 
solutions for the background evolution in the far past be- 
fore the bounce, around the bouncing point, and in the 
far future after the bounce. 



1. Solution far before the bounce 

Far before the bounce, as we assumed, (j) begins with 
a large negative value and moves slowly, while the non- 
linear term of is negligible. The energy density and 
pressure reduces to 

similarly to a single scalar field with an Ekpyrotic po- 
tential. Thus, we can choose the model parameters and 
the initial conditions of (p in order to acquire scaling so- 
lutions, namely to obtain a scale factor evolving as 

a(0~(<*-t)f ~|77,-,7|T^ , (18) 

where the subscript denotes the time where the non- 
linear term in equations (jl0p -(|13 p becomes important. 
Note that since we are considering the region where t < 
t*, we have t* — i > 0. For completeness we have also 
expressed the above solution using the conformal time 77, 
related to the cosmic time t through dt = adrj. Similarly, 
the field (/)(t) scales as 

c^{t)~--\n{t,-t),p=\. (19) 
c 

Moreover, taking the derivatives of the above expressions 
we find 



0(i)c 
H{t) 



c{U - 1) 



(20) 
(21) 



and therefore the energy density from expression 
becomes 

p(i)t<t. =3i?'cx(t-t,)-2 (22) 

Finally, using relations and p^ . for the total EoS 
using we obtain 

w ~ — — 1 = const. (23) 



2. Solution around the bounce 

Around the bouncing point Ib the nonlinear term be- 
comes important. From the Friedmann equation ([9]) we 
can express H as: 



(24) 



in which the nonlinear term becomes important and the 
last terms in pO| and pT|) can no longer be neglected. 
Without loss of generality we can keep the minus-sign 
branch in order to obtain a positive (j>{t).'^ Therefore, 
when the universe goes from the contracting (H < 0) to 
the expanding phase (H > 0), one has > y/TVoe'^'^^^ 
before and < y/2Voe'^'^^'^ after the bounce. It is there- 
fore natural to consider the solution at the bounce region 
to be 



(25) 



with a and /3 being two parameters satisfying the condi- 
tions 



a < , atB + 13^1 . 
Integrating equation ((25|) leads to 
e 2 _ e 2 = a/2Vo 



(26) 



(27) 



with to a boundary value for t, and 0o the corresponding 
value of 0. The above formula will be simplified if we set 
to in the far future when 0o becomes very large, and thus 



<='»0 



0. In this case the expression for 0(t) becomes 
0--^ln{-v^(io-i)[|(t-to) + ato + /5]} , (28) 



which leads to 



4{at + 13) 



c(to-t)[a(t-io) + 2(aio+/3)] 
2 
c 



(29) 
(30) 



It would be useful if we could approximate the above 
expressions around the bouncing point ts, namely t — 
ts = (1 — (3) /a. In this case, and neglecting the constant 



^ Note that in general the Hubble parameter could transit from 
the minus-sign branch to the plus-sign branch either before (for 
<f>initial < 0) or after the bounce (for (pinitial > 0). Thus, above 
we assume that if such transition exist it will take place after the 
bounce. However, as we will shortly see below, at late times and 
under the curvaton backreaction, expression II24I I is not exactly 
valid any more. Hence, we do not examine in detail the relation 
between the two branches. 
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terms in order to extract the pure scaling behavior, we 
obtain 



c[(/3 + ato)^-l]' 

c[{l3 + atoy--lf 
8a^ [3{l3 + ato)^ + l] jt-tB) 
c{(3 + ata- 1)^(13 + ato + l)^' 



(31) 



that is (j> exhibits a hnear behavior around ts ■ Fohowing 
the same way, one could also insert these expressions into 
equation and (|12p to straightforwardly obtain the 
approximate solutions for H{t) and w{t), respectively. 




3. Solution after the bounce 

Far after the bounce, when t approaches Iq, the solu- 
tions are still given by (P5)) - ([5n| . Thus, approximating 
them at < — > ioi and neglecting the constant terms in 
order to extract the pure scaling behavior, we acquire 



FIG. 1: The evolution of the Galileon field (j) with respect to 
t. We choose the initial conditions to he (f>i — — \/3(ln20)/3 
and 4>i = \/3/60, and the parameters to be c — 2\/3, g — 1, 
and Vo — 1/12, respectively. 



-\n{to~t) , 

c 

: 1 



cto-t ' 
2 1 
c (to - tf 



(32) 
(33) 
(34) 



In this case, inserting these expressions into ([M)) . we 
can extract a simple approximate expression for H(t) too, 
namely 



H 



1 



(io - tf 



(35) 



which leads to a Big-Rip singularity when t approaches 
to. This can be easily explained since when the non- 
linear terms become very important the last terms in 
(|TOl) and (fTTj) become dominant. Thus, when the en- 
ergy density is dominated by the 3gH(j)^ term, namely 
= p ~ 3gH(p^, we straightforwardly find that the 
Hubble parameter H will be proportional to (p^^, and ([35|) 
is verified. Finally, the total EoS can also be obtained by 
inserting these expressions into (|12p . 



4- Numerical verification 

We close the background investigation by performing 
an exact numerical elaboration in order to verify the 
above approximate expressions in the various regimes. In 
particular, we numerically solve equations © and ([T^. 
imposing as our initial conditions. In Figures [TJ [51 
[3]and[31 we respectively present (j){t), 0(t), the Hubble 
parameter H{t) as well as the EoS w{t). As we observe. 



dt^/dt 




FIG. 2: The evolution of the speed of the Galileon field 
with respect to t. We choose the initial conditions to be cj>i ; 
— V^(ln20)/3 and (pi = \/3/60, and the parameters to be c - 
2^/3, g = l, and Vo = 1/12, r 



both (pjt) and (p{t) are monotonically increasing. More- 
over, before the bounce the universe contracts with a 
scaling solution with EoS w being constant and larger 
than unity (in this specific example w ~ 3), and thus the 
scenario is free from the anisotropy problem discussed in 
section [Hi 

As time passes the nonlinear term becomes important 
and triggers the bounce, which forces H{t) to change 
from negative to positive. Note that the numerical re- 
sults confirm that (pit) has a positive value during the 
bouncing period, and thus it justifies our choice of the 
minus sign in (j24[) . After the bounce, the nonlinear ki- 
netic term makes the scalar-field energy density increas- 
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FRW metric in an Arnowitt-Dcscr-Misner (ADM) form 



N^dr + h,j{dx' + N'dt){dx^ + N^dt) , (36) 



where N is the lapse function, TV* is the shift vector, and 
hij is the induced 3-metric. One can then perturb these 
functions as: 

N = l + A, TV, = , % = a^{t)e^^5^j , (37) 

where A, ^ and C are the scalar metric perturbations. As 
usual it is useful to define the (gauge-invariant) comoving 
curvature perturbation through 



n = c 



H . 



(38) 



FIG. 3: The evolution of the Hubble parameter H with respect 
to t. We choose the initial conditions to be (j>i = — -\/3(ln 20)/3 
and (f>i = y/S/60, and the parameters to be c — 2\/3, g — 1, 
and Vo = 1/12, respectively. 




FIG. 4: The evolution of the EoS w with respect to t. We 

choose the initial conditions to be (pi = — \/3(ln20)/3 and 

4>i ~ y/3/60, and the parameters to be c = 2\/3, g = 1, and 
Vo = 1/12, respectively. 



ing, leading the universe to a Big-Rip. One can also see 
that the total EoS ■w{t) indeed indicates the bounce fol- 
lowed by the Big-Rip, verifying the analytical results that 
has been obtained in preceding paragraphs. 



and hence in the uniform (/) gauge we acquire (50 = and 
C = 7^. 

Under the above perturbation scheme, the action (|5]) 
perturbed up to second order becomes: 



5(2) = / drtd'xa'^ 



12 



(39) 



where 77 is the conformal time. In the above expression 
we have introduced the sound-speed squared c^, and the 
quantity Q-jz related to instabilities, which in our specific 
scenario read as 



i9{ 



2H^) - ]^g^^] , (40) 



,2 j,4 



(41) 



Note that according to the definition of V in p4| we ob- 



=.1.5 




B. Perturbations 

In the previous subsection we analyzed the background 
evolution of the Galileon bounce. Thus, we can now pro- 
ceed to the investigation of the perturbations, focusing 
on their stabilities. It proves convenient to foliate the 



FIG. 5: The evolution of the squared speed of sound 
with respect to t. We choose the initial conditions to be 
4>i = — \/3(ln20)/3 and = ^/3/60, and the parameters to 
be c — 2\/3, g = 1, and Vo = 1/12, respectively. 



tain 2Q-JI — T). The avoidance of ghost as well as gradient 
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instabilities requires both Q-ji and to be larger than 0. 
Using the approximate expressions of the previous sub- 
section, it can be easily verified that these conditions are 
fulfilled in our scenario. Furthermore, in order to show 
that this is the case for the exact behavior too, in Figures 
[5] and [6] we respectively depict and Qu, arising from 
the numerical elaboration of the full system. 




FIG. 6: The evolution of the instability-related quantity Q-jz 
with respect to t. We choose the initial conditions to be (j>i = 
— \/3(ln20)/3 and (pi — ^/3/60, and the parameters to be c = 
2v^, = 1, and Vo — 1/12, respectively. 



From the above analysis we can see that our model is 
stable under both ghost and gradient instabilities. How- 
ever, it generally generates a (deep) blue tilted power 
spectrum, which cannot be consistent with observations. 
Note that in this was shown for w = 1/3, there- 
fore in our present case where w > 1 the blue tilt of 
the power spectrum is even stronger. This implies that 
the anisotropy-free requirement w > 1 and the scale- 
invariant perturbation generation cannot be obtained si- 
multaneously in the current case, as was already men- 
tioned in the Introduction. For this reason, we should 
introduce an additional mechanism in order to be able 
to generate the nearly invariant perturbation spectrum, 
without spoiling the solution to the anisotropy problem. 
This can be performed by the curvaton mechanism, as 
we will present in the next section. 



IV. THE CURVATON MECHANISM AND THE 
SCALE-INVARIANT SPECTRUM 



In the previous section we showed that the Galilcon 
scenario under an Ekpyrotic-likc potential can exhibit a 
bouncing solution naturally. In the contracting phase the 
corresponding total EoS of the universe satisfies w > 1 in 
order for the evolution to be free of the anisotropy prob- 
lem discussed in section |lTl However, under this require- 
ment the corresponding perturbations, although free of 



ghost and gradient instabilities, cannot give rise to a 
nearly scale-invariant power spectrum, as it is required 
by observational data pl| . 

In [2^ it was shown that this problem can be solved by 
introducing a curvaton field cr coupled to the Galilcon (/). 
As it is usual for curvaton fields [33[ , a does not affect the 
bouncing background behavior, but it can lead to a power 
spectrum in agreement with observations. In particular, 
through a specific Galileon-curvaton coupling, the curva- 
ton field lies effectively in a "fake" de-Sitter expansion or 
matter-like contraction, and thus it can generate a nearly 
scale-invariant power spectrum. This is the idea behind 
the "conformal" mechanism [2]| and in this section we 
investigate the necessary form of the coupling functions. 

Let us consider the curvaton action as 



d^xV^j[-Tmdaf - 2g{4>)W{a)] , (42) 



allowing for the most general form of coupling between a 
and 0. The functions J^{(j)) and g{(t>) depend on 0, while 
W{a) is the potential for a. Variation of action with 
respect to a gives its background evolution equation as: 



-cro 



Wa„ = , 



(43) 



where uo is the background value of a and W^q corre- 
sponds to dW/d(j\cro- In the following, and up to the end 
of this section, we omit the subscript "0", denoting the 
background by a simple cr, unless explicitly mentioned. 
Additionally, the energy density and pressure of a can be 
respectively written as: 



Pa = \Fa^ +gw , P^^ ^Ta^ 



gw 



(44) 



We now proceed to the examination of the perturba- 
tions generated from the curvaton field. Perturbing it by 
da and defining u = a\f^5a, we can extract the corre- 
sponding perturbation equation as a second-order differ- 
ential equation: 



F z 







(45) 



where the prime denotes derivative with respect to the 
conformal time 77, Waa = d'^W {a) / {daY is the second 
derivative of the potential with respect to a, and wc have 
defined z = a\lT . As usual, the power spectrum gener- 
ated by (5(7 is defined as: 



(46) 



Therefore, we deduce that the condition for obtaining a 
scale-invariant power spectrum of ba is: 



a^gwg 

T 



z 
z 



l»7* - '71' 



(47) 



There are several ways to satisfy the condition (|47p . 
The simplest one is to set W[a) = 0, in which the first 
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term in the above condition disappears, and thus we need 
just to suitably choose T in order to obtain z" j z ~ 2\r]^ — 
r]\~'^. Alternatively we can incorporate the effects of both 
the kinetic and the potential terms of ct. In the following 
subsections we consider these cases separately. 



A. W{a) = 

Under W{(j) = 0, the condition (jT7|) of obtaining scale- 
invariant power spectrum becomes 

z cx - ?7|^ or \t]^ - r/p-^ . (48) 

For the case of z cx |r/, — r^p wc obtain 

(5a = ^~fci, k-^Tj,-Tj\-\ (49) 

the latter of which dominates over the former. Therefore, 
using expression (|46p we can obtain the power spectrum 
as 



(50) 



As we observe, the spectrum is indeed scale-invariant but 
it has an increasing amplitude. 



On the other hand, for the case of z oc 



acquire 



Sa ^ - ^ k 2 ^ fc2 177* — 77I 



ri\ we 



(51) 



the former of which dominates over the latter. Therefore, 
using (|46p wc can obtain the power spectrum as 



(52) 



In this case the spectrum is scale-invariant and moreover 
it is conserved on super-horizon scales. 

The absence of W{a) leads to an absence of G{<j)) too. 
Thus, we only need to suitably determine the form of 
J^{4>) according to the condition (|^5)) . Note that far be- 
fore the bounce we have already assumed the scale-factor 
ansatz (fT8|) . and therefore ([48]) requires just 



2(2-.3p) 2 

CX jry* — 77I i-p or — ?]\p-^ 



(53) 



Furthermore, from the ansatz solution (|19p for (j) and the 
above expressions we deduce that the suitable choice of 
the form of in terms of (j) might be 



cx e2(3-c^)0/c 



or e 



c<p 



(54) 



where we have made use of the relation p = (note 
that since in contracting phase w 3> 1, we have < p <C 

We close this subsection by examining the backreac- 
tion of a field on the background evolution, since al- 
though the curvaton is necessary for the correct per- 
turbation generation we would not desire it to spoil the 



background bouncing behavior itself. In the contraction 
region where t < t^, the background energy density of 
the system scales as pt<t, ^ {t* ~ t)~'^, as it was found 
in (P^ . On the other hand, the evolution equation of 
the curvaton field ([43]) gives & ^ a^^T^^, and thus its 
energy density in (|^^ becomes — J'cP' jl ^ a~^J-~^ . 
Since we know the time-dependence of the scale factor 
from and the time-dependence of J- from ([55]) . we 
straightforwardly deduce that for the solution branch 
where z cx \r\^ — we obtain ^ (t* — i)^*, while 
for the solution branch where z cx |?7* — ?7|^^ we ac- 
quire p„ - (t* - t)2(i-3p) ^ _ if (in the last step 
we used that p ^ 1). Therefore, our analysis indicates 
that in the solution branch where z cx \t\* ~ (with 
T cx e^'y^~'^ )<i>l'^"^^ one has to suitably tune the initial 
conditions in order for the curvaton not to destroy the 
background bouncing behavior. However, in the second 
solution branch where z cx jjy* — 'r]\~^ (with T cx e'^'^), the 
energy density of the curvaton field grows slower than 
that of the background, and thus the background bounc- 
ing evolution is not altered by the backrcaction of the 
curvaton. 



B. W(o-) / 

We now examine the case where the curvaton poten- 
tial is non-zero, in order to investigate its effect on the 
perturbation generation. Without loss of generality and 
for simplicity we assume that J- is approximately a con- 
stant (we set T = 1), although extension to general T is 
straightforward. 

In order to see what condition (|T7|) gives in this case, 
we recall that at the early stage of the bouncing phase, 
where the perturbation 8(j is generated, the scale factor 
evolves according to ^5]) . and therefore since z = a\fT 
we obtain 



z" a" Pip 



- 1 



1 



z a 1-p \ l-p J \T]-r]^f ' 
Furthermore, introducing — a(r]^) we can write 



(55) 



(56) 



Inserting equations (|55|) and ((56)) into (|47)) we deduce 
that the condition for obtaining a scale-invariant power 
spectrum of 5a becomes 



3p-2 
(1^ 



\v*~v\ 



(57) 



As a specific example we consider the well-studied case 
of a quadratic potential, namely W{a) ~ m^(T^/2, in 
which case W^a = const. Hence, condition (|57p . using 
also the (/)-evolution from (|T9)) . gives 



G oc e 



C(p 



(58) 
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Therefore, we extract that the field perturbation 5a 
scales as: 



0(7 = — 



3p-2 



(59) 



We mention that since we assume J-' = 1 (that is z = a) 
the "fake" effect is absent, and thus the dominating mode 
of the perturbations is always the growing mode. 

We close this subsection by examining the backreaction 
of a field on the background evolution, since we would 
not want the curvaton to spoil the background bouncing 
behavior itself. The cr-evolution equation becomes 



a + 3H& + mlGa = 



(60) 



Since according to relations ([58|) and ([T9| we have G 



(t* — t) ^, we can write G = Go{t* — t) with 



Go being an arbitrary constant. Therefore, equation ((60 
accepts the solution 



(t*-t)5[l-3P±V(l-3p)^-4- 



iGo] 



(61) 



Finally, substituting it into (|44p for the curvaton energy 
density we obtain 



{t^ - i)[l-3p±V(l-3p)2-4m2 501-2 



(62) 



Comparing the background energy-density evolution 
with the curvaton energy-density evolution (|62p we 
deduce that the requirement for the latter to grow slower 
than the former is to have 1 — 3p±-\/(l — 3p)^ — Am^Go > 
0. However, since in order to solve the anisotropy prob- 
lem we focus on w > 1 (or equivalently p < 1/3), then 
provided Go > ^ the above requirement is always satis- 
fied. Therefore, in the scenario at hand the energy den- 
sity of the curvaton field will never dominate over the 
background evolution, that is the background bouncing 
behavior will not be destroyed by its backreaction. 



V. RECONSTRUCTING THE EXACT 
SOLUTION AROUND THE BOUNCE 

In the previous sections we constructed the Galileon 
bounce scenario free of the anisotropy problem, in which 
we added the curvaton field in order to obtain a scale- 
invariant power spectrum of primordial perturbations 
generated in the contracting phase. Additionally, we 
showed that under soft requirements on the choice of 
T{(f>) or Giff)), the backreaction of the curvaton field will 
not alter the background bouncing evolution. 

However, after the bounce the effect of the curvaton 
field can be significant, and in particular it can regu- 
larize the universe evolution in order not to result to a 
Big-Rip. In order to examine what classes of coupling 
functions can provide this overall behavior, in this sec- 
tion we semi-analytically reconstruct them following the 
"inverse" procedure (30| , in which we impose as input the 
desired bouncing scale factor, reconstructing suitably the 



various function in order to correspond to a consistent 
and exact solution of the full system of equations. 

The complete action of the Galileon-curvaton system, 
consisted of both ([5]) and (|^^ . can be written as: 



5, 



total 



1 



1 



+ |V^0V'^0n0- J-(0)(9a)2 -2g(0)VF(a)J .(63) 

Note that although matter and radiation could be in- 
cluded straightforwardly, in the above action we have 
neglected them in order to examine the pure effects of 
the Galileon-curvaton evolution. Thus, the cosmologi- 
cal equations in the FRW metric arc the first Friedmann 
equation: 

and the evolution equations for the two fields, namely 



5J-(0) 



90 J-(0) 



g(0) dW{a) 
J"(0) da 



(65) 



and 



1 + QgHcj) + -g^ 



2 06 ^ 



dG{<l>) 



3 
2 

dV{(t)) 



2H + gcf) 
, 



F{(t))a'^ - ?>gH4r 



W{a) 



(66) 



respectively. 

One could solve the above equations fully numerically, 
imposing specific ansantzes and initial conditions, how- 
ever doing so he does not have control of what functions- 
ansantzes, parameter choices and initial conditions, lead 
to bouncing solutions. That is why the aforementioned, 
semi-analytical, "inverse" procedure, where the scale fac- 
tor is imposed a priori^ is better and more appropriate 
for the analysis of this section, allowing for a system- 
atic control on the conditions of the bounce realization. 
We mention that since there are more unknown functions 
than equations, one can in general always reconstruct the 
desired evolution. 

Let us impose a desired bouncing scale factor a{t) as 
an input, by which H{t) is also known. Furthermore, 
we consider V{(j)) as usual, and we impose (j){t) at will 
too. Thus, only three free functions remain, namely 
J-{4)), Gifp) and W{a) which must be derived by the equa- 
tions, along with the solution for a(t). Since there are 
three independent cosmological equations, namely equa- 
tions ([M)) . (p5|) and we must also impose by hand 
one more of the above four functions. We prefer to set 
W{a) = 0, since this was one (simpler) case that was 
approximately analyzed in subsection IIV Al (however one 
could easily consider other W{a) forms too). Such a 
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choice simplifies things since the function G{<t>) also dis- 
appears from the equations, and therefore the cosmo- 
logical equations (|64 |) -(|66 |) are considered as differential 
equations for cr(i) and -Fit). Thus, after obtaining the so- 
lution, and since we know 4'{t), we can reconstruct J^{4>)- 
Equation ([64|) can be algebraically solved in order to 
obtain as 



1 



m{tf~2vm))^^\t) 

^QgH{t)^\t)\ . (67) 



Substituting this into (|65p gives a simple first order dif- 
ferential equation for F{t) of the form 



h{t,T{t),i-{t)) = 



(68) 



which can be easily solved. Thus, from the solution of 
J^{t) and the known (j){t) we can reconstruct J-{(j>)- 

We mention here that the above procedure holds for 
every input functions, with the only requirement being 
the obtained from relation (p7)) to be positive, other- 
wise there is no solution that can correspond to the input 
functions. With the above semi- analytical procedure one 
has full control on how to choose the model parameters 
and the initial conditions in order to get a positive in 
([67| . On the other hand, if one tries to solve fully nu- 
merically the three equations simultaneously, it 
is very hard to determine the model parameters and the 
initial conditions in order to get a consistent bouncing 
solution. 

In order to apply explicitly the above reconstructing 
procedure, without loss of generality we choose a bounc- 
ing scale factor of the form 



a(t) = as 



1 3 . 

1 + -ur 



1/3 



(69) 



where is the scale factor at the bouncing point and uj is 
a positive parameter which describes how fast the bounce 
takes place. The above ansatz presents the bouncing 
behavior, where t varies in the bounce region, that is 
between the times ts- and tB+-, with t = ts ~ the 
bouncing point, however one could use at will any other 
bouncing ansatz. Straightforwardly we find 



Hit) 



(1 -h3wtV2) 



(70) 



and thus the universe is free of a Big-Rip after the 
bounce. 

For the field and the potential V^((/>), enlightened 
by the analysis of subsection IIII Al and without loss of 
generality, respectively we assume 



and 



0(i) = \n{t - ti) 



(71) 



(72) 




FIG. 7: (Colored online) The solution for the coupling func- 
tion T{t) and the imposed Galilean field (t){t), under the im- 
posed bouncing ansatz S69\) . We choose the parameters as 
g = 1, Vo = 1/12, c = 2V^, cj = 1, (f)i = 0.1, ti = -100, and 
tB± = ±1. 




0.456 0.457 0.458 0.459 0.460 0.461 0.462 0.463 0.464 0.465 



FIG. 8: (Colored online) The reconstructed coupling function 
J-{ij>) under the imposed bouncing ansatz \69\) . using Fig. 
We choose the parameters as g — 1, Vo — 1/12, c = 2v3, 
= 1, 0/ = 0.1, ti = -100, and tB± = ±1. 



while as we mentioned we set W{a) = 0. 

We follow the procedure described above, and for the 
model parameters we choose g = 1, Vq = 1/12, c = 2^3, 
a; = 1, 0j = O.I, ti = -100, and tB± = ±1. In Fig. [7] 
we depict the solution for F{t) and also the known 0(t) 
from (|7ip . and in Fig. |S] we present the corresponding 
reconstructed F {(/)). Finally, for completeness in Fig. [9] 
we show the solution for (T(t). 

We close this section by mentioning that in principle 
one could think of other reconstructing procedures, for 
instance setting F{(j)) and reconstruct G{<P) and W{a), 
or even setting (j){t) and a{t) and reconstruct F{(j)), G{<l>) 
and W{(j). So there can actually be many possibilities 
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FIG. 9: (Colored online) The solution for the curvaton field 
a{t), under the imposed bouncing ansatz S69}) . We choose the 
parameters as g = 1, Vo = 1/12, c = 2\/3, uj — I, (fii = 0.1, 
ti = —100, and tB± = ±1- 

to realize the non-singular bounce. 

VI. CONCLUSIONS 

Bounce cosmology is an interesting paradigm since it 
alleviates the Big-Bang singularity problem. Addition- 
ally, it can solve the Big-Bang problems, and nearly scale- 
invariant primordial perturbations can be incorporated 
too. These features make bouncing cosmologies success- 
ful alternatives to inflation. However, there are many 
detailed issues that should be carefully addressed dur- 
ing the establishment of bouncing cosmology, and in the 
present work we tried to confront some of them. 

First of all, the NEC violation, which is required for 
the bounce realization, may bring ghost degrees of free- 
dom. In the above analysis we were based on the Galileon 
scenario, which is a higher-derivative construction free of 
ghosts, and thus we obtained a bouncing evolution free 
of ghost and gradient instabilities. 

However, there is a secondary problem that may dis- 
turb the bounce construction, namely that in the con- 
tracting phase even a tiny anisotropic fluctuation from 
the totally isotropic FRW geometry will be radically en- 
hanced and destroy completely the FRW evolution. The 
solution of this "anisotropy problem" requires the total 
EoS of the universe to lie in the regime w > 1 in the 
contracting phase. Thus, starting from [2^ where the 
anisotropy problem was present, in this work we were 
able to solve it and obtain w > 1 by considering an 
Ekpyrotic-like potential with negative value. This is one 
of the main contributions of the present article. 

The above solution of the anisotropy problem through 
a large EoS has an undesired effect, namely it spoils 
the generation of a nearly scale-invariant power spec- 
trum, that a scalar with smaller EoS can bring through 



adiabatic perturbation. Therefore, in order to still be 
able to produce a power spectrum in agreement with 
observations, we additionally introduced in the scenario 
a second, curvaton field, coupled to the Galileon one, 
which can indeed generate the desired perturbations in 
an isocurvature way. In our analysis we presented this 
mechanism in general, and we analyzed explicitly two 
specific examples where nearly scale-invariant perturba- 
tions are generated. Finally, we examined the conditions 
under which the curvaton field does not cause a signifi- 
cant backreaction on the background bouncing behavior 
caused by the Galileon field. 

Furthermore, the curvaton field, apart from the gener- 
ation of the desired perturbations, has another important 
role, namely after the bounce it can regularize the back- 
ground evolution in order to avoid a Big-Rip singularity, 
which is caused by the Galileon field itself. In partic- 
ular, although the curvaton backreaction is not signifi- 
cant at the background level before the bounce, during 
and after the bounce it becomes important and changes 
the background evolution. In order to see this effect we 
performed a semi-analytically "inverse" analysis, recon- 
structing suitably the desired bouncing evolution of the 
scale factor, which is free of a Big-Rip without any fine- 
tuning. 

We close this work by mentioning that there could be 
other possibilities to avoid the Big-Rip singularity. For 
instance, an alternative evolution after the bounce would 
be to assume that the Galileon and curvaton fields de- 
cay to standard model particles [13]. In this case the 
decaying Galileon energy density cannot trigger the Big 
Rip anymore, and additionally it can produce the matter 
content of the universe. Such a detailed analysis of the 
post-bounce evolution, and its relation to the subsequent 
thermal history of the universe, lies beyond the scope of 
the present work and it is left for a future investigation. 

Note added: After completing our manuscript, we 
came to know that studies of bounce cosmology aiming 
to the same issue has been done in [3^, in which simi- 
lar results are obtained for different (conformal) Galileon 
models but with the same (Ekpyrotic-like) potential. 
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